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We propose an effective Ginzburg-Landau energy func- 
tional for superfluid helium that is motivated by the ideal 
Bose gas model. The parameters of this phenomenological 
model are adjusted to known properties of the bulk system. 
The motivation and the field of applications of this model 
are quite similar to that of the 5' theory by Ginzburg and 
Sobyanin. The correlation length, the depression of the tran- 
sition temperature in helium films and the density profile at 
a boundary are calculated. 



I. INTRODUCTION 

In the Ginzburg-Landau (GL) model for liquid "'He the 
superfluid density gs is identified with the expectation 
value of the square of the complex order parameter field 
*(r), 



(1) 



Here m is the mass of a helium atom. In the vicinity of 
the A-transition, the energy functional is expanded into 
powers of the field ^ and its derivatives and of the rela- 
tive temperature t = T /T\ — 1. In its simplest form the 
GL energy functional takes into account the terms with 
|V*p, t • and l^'l''. 

The GL model for superfluid helium has been intro- 
duced in 1958 by Ginzburg and PitaevskiiEl. It leads to 
critical properties (like Qs oc \t\) that do not agree with 
the experiment. This failure is due to the fact that in 
this model the fluctuations diverge (relative to the mean 
value) when the transition point is approached. 

The theoretical remedy for this break-down of the GL 
model is well-known: The GL ansatz may be considered 
as a sensible approximation for flnite volumes of the sam- 
ple only. Using scale transformations to larger and larger 
volumes renormalizes the coefficients in the energy func- 
tional. The renormalization group theory (RGT) leads to 
nontrivial critical exponents that are no longer in vari- 
ance with experiment. RGT is considered as the best 
present-day theory which allows in principle the calcula- 
tion of all relevant quantities. 

For practical reasons it is, however, useful to have a 
phenomenological GL like model that is in accordance 
with the bulk properties and in particular with the ex- 
perimental critical exponents of superfluid helium. Such 
an effective GL model allows to address a large variety 
of problems in a simple and straight-forward way. 



An effective GL model of this kind has been put 
forward Wider the name theory" by Ginzburg and 
Sobyanincl^Q. The references ||-^ are comprehensive pa- 
pers addressing the foundation and applications of such 
an approach. For the historical development of the field 
we refer to the corresponding discussion and the bibliog- 
raphy in these papers. 

The ^ theory by Ginzburg and Sobyanin (abbreviated 
by GS in the following) is based on Eq. (|l|) , too. The GS 
energy functional contains the terms jV^'p, • l^'P, 

|i|2/3.|vj,|4 g^j^^ |^j>|6^ The temperature dependences of the 
coefficients are chosen by hand such that they yield ap- 
proximately the right critical exponents (like Qs oc 
for example). In so far as the actual critical exponent of 
Qs deviates from = 1/3 this model is an approximation. 
For the sensibility and usefulness of such an effective GL 
model and its relation to RGT we refer to the thorough 
discussion by GS0 

It is generally accepted that there is an intimate 
connectionOij between the Bose- Einstein condensation of 
the ideal Bose gas (IBG) and the A transition in Hquid 
helium. In Ref. || we have argued that there are good 
reasons to uphold the IBG critical exponent j3 — 1/2 
of the condensate density qq = m(|^'p) ~ |t|. This 
implies qq < Qs just below the lambda point. Conse- 
quently, we have proposed that the superffuid density is 
composed by the condensate density and a coherently 
comoving density ^?coh (corresponding to lowlying non- 
condensed states): 



Qs — Qq + Qcoh 



where 



£'0 



(2) 



The arguments in favor of this composition of Qs and the 
resulting expression for Qs are reviewed iiH Sec. II A. The 
expression for Qs provides excellent fitsu to the experi- 
mental temperature dependence of the superfluid density. 

Based on the starting point ^ instead of Eq. (|l|) we 
construct in this paper an alternative effective GL model. 
Because of go = m ~ \t\ we may use the standard 

form for the Landau part of the energy functional (i. e. 
the terms t • j^fp and |5'|'*). The kinetic energy contribu- 
tion (the term with jV^'p) gets, however, an additional 
factor Qs/ Qo due to the coherent comotion. Effectively, 
this factor Qs/qo ^ 1^1"^^'^ damps the critical ffuctuation 
such that our GL ansatz becomes scaling invariant. 

From a pragmatic point of view our effective GL model 
(with the terms • jV^-p, t ■ \^\'^ and |*|4) is quite 

comparable to the 5" theory by GS (with the terms |V^p, 
|^|4/3 . |^|2^ |^|2/3 . |^|4 |^|6)^ j,^^ ^^^^^ difference is 

that in our model the Landau part is of its standard form 
(in accordance with the IBG), and that specific reasons 
are given for the occurrence of a fractional power of \t\. 
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This paper is organized as follows. Section |l| presents 
the reasons for the composition (H) of the superfluid den- 
sity, the actual expression for Qs, and the parameters ob- 
tained by a fit to the experimental temperature depen- 
dence of Qs- Section [II displays the energy functional of 



our effective GL model, and discusses its basic features. 
Section IV investigates the temperature dependence of 
the correlation length and the size of the critical fluctu- 
ations. The depression of the transition temperature in 
helium films is calculated in Sec. 0. Section fvl derives 



the spatial variation of the superfluid and the condensate 
density at a boundary. 



II. COMPOSITION OF THE SUPERFLUID 
DENSITY 

A. Motivation 

The comparison between liquid "^He and ^He on one 
side and the ideal Fermi gas and the ideal Bose model 
(IBG) on the other side strongly suggests that there is 
an intimate relation between the Bose-Einstein conden- 
sation (BEC) and the lambda transition in liquid ''He. 
It is, moreover, inviting to identify 5" in Eq. (|l|) with 
the IBGUeondensate wave function. This identification 
explaing3^Q a number of|-experimental findings of which 
the most important oneu is that a superfluid current is 
irrotational. The IBG and this identification are, how- 
ever, in conflict with the experimental critical exponent 
ly K 1/3 of ~ jip" (as compared to the IBG with 
/3=l/2for(|*|2)^|t|2/3). 

The standard solution of this conflict appears to be the 
renormalization group theory (RGT). The RGT yields in- 
deed an explanation of the experimental value of ~ 1/3. 
It does, however, not resolve the conflict with the IBG 
value (3 — 1/2: A renormalization is appropriate for 
the Landau value /3l = 1/2 but not for the IBG value 
(3 = 1/2. The reason is that the IBG value is obtained by 
the exact evaluation of a partition sum. This exact eval- 
uation implies a summation over arbitrarily small mo- 
menta (or, correspondingly, arbitrarily large distances). 
Therefore, the reasoning behind the renormalization pro- 
cedure (analytic Ginzburg-Landau ansatz for finite re- 
gions or blocks, and subsequent transformation to larger 
and larger blocks) cannot be applied to the IBG energy 
functional. 

Within the framework of the ideal Bose gas model 
we have proposedB to resolve the conflict /3 = 1/2 and 
v ~ 1/3 by the assumption that noncondensed particles 
move coherently with the condensate. This means that 
we no longer identify the condensate with the superfluid 
fraction; the condensate is only part of the superfluid 
phase. The coherent comotion is constructed in the fol- 
lowing way: The condensate wave function may be writ- 
ten as 



where (pQ{r) and $(r) are real functions. The velocity of a 
potential superfluid flow is given by Vg = (?i/m)V$. The 
single particle function </3k(r) of the noncondensed states 
can be chosen as real functions. If now the lowlying states 
(up to a coherence limit fccoh that is to be determined 
later) adopt the same phase factor. 



ipk{r) ^k(r)exp[i$(r)] , (fc < fccoh) 



(4) 



then the corresponding particles contribute to the super- 
fluid current density. Consequently, the superfluid den- 
sity is made up by the condensate density go and the 
density Pcoh of the coherently comoving particles: 




go + gcoh 
Q 



(5) 



Here (n^) are average occupation numbers of the (slightly 
modified) IBG, and the prime at the summation symbol 
excludes the condensed particles from this sum. 

We emphasize that the construction (^ to (P) pre- 
serves the basic properties of the superfluid flow (like its 
irrotational behavior). In section HI C| we discuss the 
physical meaning of fccoh and present an argument why 
the lowest single particle levels might adopt the same 
phase factor as the condensate. 



B. Functional form of the superfluid density 

We determine the asymptotic form of the expression 
(^). The IBG expectation values {uk) are given by 



1 



1 



{rik) 



exp [(efc - ^l)/kBT] - 1 exp(K2 +t^)-1 



(6) 

Here /i is the chemical potential, Ek = h^k^ /2m are the 
single-particle energies, and fee is Boltzmann's constant. 
In the last expression we introduced the dimensionless 
momentum 



Alkl 



where A 



An 



\/2t: rnksT 



(7) 



The chemical potential /i has been expressed by the di- 
mensionless quantity = — /i/fceT in Eq. (^. We 
expand this r into powers of the relative temperature 
t = T/Tx - 1: 



Tit) 



-/i _ ( at + bt^ + . 
feT " 1 a'\t\+b't^ 



(t>0) 
{t<0) 



(8) 



l'(r) = v3o(r)exp[i$(r)] 



(3) 



For t > the particle number condition N = X]k("'k) 
determines the temperature dependence of T{t) and in 
particular the coefficients a, 6,. . . , for example a — 
3C(3/2)/(47ri/2)^ (jjgj.g ^ Er"""" denotes Rie- 

mann's zeta function). For t < the IBG yields r = 0. 



2 



As a modification of the IBG we admit nonvanishing co- 
efficients a', • By a' 7^ we introduce a phenomeno- 
logical gap between the condensate level and the noncon- 
densed levels. This modification does not affect the BEG 
as the most important feature of the IBG. It is necessary 
in order to reproduce the experimental temperature de- 
pendence of the superfluid density by Eq. (^. In view 
of the successful roton picture it is not surprising that a 
gap is required for a quantitative description. 

For t < the particle number condition N = (no) -t- 
Y^'ki^k) leads to 



go 
Q 



N 



f\t\ + 0{t') 



where 



/ 



3 2^ a' 

2^ am 



(9) 



(10) 



An evaluation of the coherently comoving density in 
(|) yields 



Scoh 



1 

N 



4:Kr 



higher order . (11) 



V^C(3/2) 

The leading experimental behavior is approximately 



El 
Q 
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(12) 



Deviations from the exponent 2/3 are neglected in our 
phenomenological model. Equation (021) corresponds to 



where 



ao 



4ko 



0FC(3/2) 



(13) 



(14) 







ft = PO + Pcoh 
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FIG. 1. Decomposition of the superfluid density ps into 
the model condensate density po and the coherently comov- 
ing density pcoh as a function of the temperature T. For 
T ^ T\ the comoving density pcoh oc \tf^''^ is the dominant 
contribution to the superfluid density. 

The model condensate density Qq may be multiplied 
by a depletion factor (roughly 0.1) in order to become 
comparable to the measurable condensate density. The 
temperature dependence of the resulting condensate den- 
sitji-jCompares well with experimental results by Snow et 

ai.y. 

In this paper we need only the coefhcients of the lead- 
ing terms in Eqs. (^ and (|l2|). The fit of the theoretical 
expression to the^xperimental data yields the following 
parameter valuesa 



a' = 3.03 , Ko = 2.70 



(15) 



According to Eqs. ( [Tol) and (|l4| ) these parameters are 
equivalent to 



/ = 5.60 , ao = 2.33 



(16) 



These two parameters are used as an input for our energy 
functional. 



C. Results 

An excellent agreement with the experimental tempeij^ 
ature dependence of the superfluid density is obtained^ 
if higher order terms in Eqs. (|^) and (|ll|) are included. 
The resulting fit formula contains four parameters as do 
comparable standard fits. 

This fit defines the temperature dependence of the con- 
tributions Pcoh a-nd gnto the superfluid density. The re- 
sult is shown in Fig. |l| that has been taken from Ref. 0. 
For |i| < 0.01 one finds Qs ~ Qcoh (because of £<coh ^ \t\^'^ 
and Qq ^ \t\). Due to the gap between the condensed and 
the noncondensed particles, the contribution pcoh is sub- 
stantial only in a narrow region (|i| < 0.2). For \t\ > 0.2 
one obtains Qs ~ go- 



III. EFFECTIVE GINZBURG-LANDAU 
FUNCTIONAL 

A. Definition 

The equilibrium value of Gibbs free energy G depends 
on the temperature T, the pressure P and the total par- 
ticle number N. In a GL like model, the thermodynamic 
potential G depends, moreover, on the order parameter 
field The minimum of G with respect to ^ determines 
its expectation value. The equilibrium value of the po- 
tential is then G(r, P, N) = (G(*)). 

The number of particles is conserved in the considered 
applications. The pressure dependence will not be dis- 
cussed explicitly; it enters via Tx(P) and via the pressure 



3 



dependence of the parameters (16) that may be deter- 
mined by the fit to the experimental superfluid density 
Qs{T,P). The foUowing considerations concentrate on 
the t and dependences. 

Equations (|2) and (|^) yield the bulk equilibrium value 



w(|«'(bulk)p) 



N 



I\t\ 



(17) 



Here v — V/N is the volume per particle. The order 
parameter field 4'(r) will, in general, depend on the co- 
ordinates. The thermodynamic potential is, therefore, 
written as an energy functional. 

We split the thermodynamic potential into several 
parts: 



G — Gq + Gql — Go + Gfluct + Gl 



(18) 



The contribution that is not directly related to ilie order 
parameter field is denoted by Go- Following GSHu we as- 
sume that the approximately logarithmic singularity of 
the specific heat is due to this term. Our investigation 
is restricted to the Ginzburg-Landau (GL) contribution 
Gql that is due to the order parameter field. This con- 
tribution may be further split into the Landau part Gl 
and a part Gfluct that contains derivatives of the field 'J. 

Near the transition point the energy functional is ex- 
panded into powers of the order parameter field: 



Ggl = 



(fir 



{A\\/^\ 



G|*|4 



(19) 



First we consider the terms without derivatives, i.e. the 
Landau part Gl- The most simple standard Landau 
ansatz {B = ht, C = const, and no higher order terms) 
leads to the IBG like behavior (^7|). The coefficient in Eq. 
(0) is reproduced if we set Gl oc t |*|V/-hw|'J'|V(2/2 
For the determination of the prefactor we follow GS" 
who argue that Gl should account for the jump of the 
specific heat at the lambda point. This leads to 



to be replaced by Qs times the squared velocity, 
implies 



&(bulk) ao|i|2/3 



go (bulk) 



f\t\ 



This 



(22) 



For the last expression we used Eqs. (|[) and (p^). By 
R ^ 1^1^^^'^ we admit a fractional power of the relative 
temperature in the coefficients (similarly as in the ^E" the- 
oryby GS). 

From Eqs. (|^), (^) and ^ we obtain the central 
expression for the energy functional: 



Ggl 



= Icfir 



aoA2|V*|2 Acp (tYi 



47r/|t|i/3 



/ 



m 

2/2 



(23) 



Here Ac denotes the thermal wavelength at the critical 
temperature: 



Ac 



47r ^/2mkBTx 



1.67A . 



(24) 



For the numerical value we inserted the mass m = m(He) 
of a helium atom and T\ « 2.17K. The other coefficients 
in the energy functional are given by 



ao = 2.33 , / 



.60 



Acp 



2.77 



(25) 



The first two parameters are determined by the fit of 
Eq. (||) to the experimental superfluid density (Sec. |II C] ). 
The value for Acp is taken from Table V (third line) of 
Rcf. 0. To some extent, the parameter values depend 
the other terms in the corresponding fit formula used. 

We usepreferably dimensionless quantities: Both sides 
of Eq. ( p3| ) are dimensionless. The integral J(Pr\'^\'^ 
is dimensionless and of order TV. The product AcV is 
dimensionless, too. The occurring numerical factors (ao, 
/ and Acp/fce) are of order 1. 



Gl 
kBT\ 



- d 



Acp 



/ 



~2p 



(20) 



where Acp is the jump of the specific heat per par- 
ticle. The bulk expectation value of the r.h. s. for 
i < is -ViAcp/kB)\t\yi2v) yielding (Gl/TV) = 
-AcpTx\t\^/2. For t > we have (Gl/TV) = 0. The 
second temperature derivative of these expectation val- 
ues reproduces the jump of the specific heat. 
The contribution Gfluct may be written as 



Gflu 



1 



ksTx kBTx 



2m 



d^r — R\W'i>\ 



(21) 



In a standard GL ansatz the factor R equals 1. According 
to Eqs. (^)-(^) the density Pcoh moves coherently with 
the condensate density qq. The term jV^Pp corresponds 
to the product of go times a squared velocity. This has 



B. Discussion 

1. Comparison to other approaches 
We consider the functional 

Ggl = J d^r {A |V*p + B + C |*|-* + D \^\ 



(26) 



and compare the coefficients of our model with that of 
a standard Ginzburg-Landau ansatz (GL) and of the \E' 
theory (GS): 





A 


B 


c 


D 




1 




1 





G^ 


1 


|i|'/3 


|t|2/3 


1 


this work 


\t\-"' 


t 


1 






(27) 
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The entrance 1 stands for a nonvanishing constant. For 
GS the IvPl^ term must be included since it scales with t'^ 
(as do the terms and |tp/3|4f|4). Because of this 

term the GS functional contains an additional adjustable 
parameter. In the standard GL as well as in our approach 
the l^*!^ term is a higher order term. 

Due to the connection to the IBG, our Landau part 
(coefficients B, C and D) coincides with the standard 
GL form. The occurrence of the dependence in 

the fluctuation term (coefficient A) is due to the coherent 
comotion of noncondensed particle with the condensate. 

As far as the results are concerned, our approach is 
more similar to GS than to GL. In particular the corre- 
lation length scales with as in the 5" theory, and 
not with as for GL. 

2. Choice of coefficients 

The prefactor of the Landau part (Effl) has been cho- 
sen such that the expectation value (GT) reproduces the 
jump Acp of specific heat at the lambda point. This 
choice may be criticized for two reasons: (i) The fluctu- 
ation part (Gfluct) leads to a P term for t < 0, and will 
contribute to the jump, too. (ii) There might be a con- 
tribution to the jump from the part (Go) in Eq. (|l8|). By 
assuming that the jump is solely due to the Landau part 
(Gl) we follow for simplicity GS. The alternative is to re- 
place Acp/fce by a unknown parameter. Adjusting this 
parameter later (for example to the correlation length 
amplitude) would yield a value rather close to Acp/fce. 
This might be considered as a pragmatic reasoning of the 
present choice. 

The factor R — g^/ go introduced in Eq. (El|) has been 
motivated by the coherent comotion (Eqs. ^ to (||)) of 
noncondensed particles. This comotion refers to a poten- 
tial superfluid flow with Vg = (?i/m)V$(r). The moti- 
vation for the i?-factor has thus been based on a rather 
specific coordinate dependence of \E'(r) for t < 0. It is, 
however, the spirit of a GL like approach to consider 
only one functional for the complex ^'(r) field, i. e. to use 
the same coefficients for all occurring coordinate depen- 
dences, and for t < and t > 0. Alternatively, one could 
introduce different coefficients for the coordinate varia- 
tion of the amplitude and the phase of the \E'(r) field; this 
would result in a more complicated model. 

We have used the asymptotic form R ^ in the 

functional (|2^). This is justified only as long as the lead- 
ing asymptotic expressions gs |tp/'^ and go ^ \t\ are 
good approximations. This is the case for 

|i| < 0.01 (range of applicability) . (28) 

The fit! of Eq. (||) to the experimental superfluid pro- 
vides the temperature dependence of the ratio R(t) in a 
wider range. Using this more general function R{t), our 
approach might be extended into the less asymptotic re- 
gion, for example to 0.1 > \t\ > 0.01, where a GL like 



approach is still meaningful. The emphasis of the present 
work is, however, on the temperature range (^). 

3. Coordinate dependence of the order parameter 

We discuss the coordinate dependence of the order pa- 
rameter field and its relation to the condensate and su- 
perfluid density. The order parameter field may be writ- 
ten as 

*(r) = |*(r)| exp[i$(r)] . (29) 

We are thus dealing with two real functions, the am- 
plitude |5'(r)| and the phase $(r), respectively. These 
functions are determined by the field equations that fol- 
low from the energy functional ( |23| ) and that are supple- 
mented by suitable boundary conditions. 

According to Eqs. (^ to (^) all particles contribut- 
ing to the superfluid density adopt the phase ^{r). The 
common phase factor has, however, no implication on the 
coordinate dependences of £io(r) and Pcoh(i")- 

The order parameter field is directly connected to con- 
densate density, Eq. (||). This imphes 

eo(r)-m(|vI/(r)p) . (30) 

It is not a priori clear how to determine the coordinate 
dependence of gicoh(r)- Eq. ( p^ ) fixes the bulk relation 

ft(bulk) mR{t) (|«'(bulk)p) (31) 

between the ^ field and the superfluid density, only. The 
coordinate dependence of gis(r) versus that of |^'(r)p will 
be discussed in Sec. ^ for the important case of the 
behavior at a boundary. 

C. Fluctuations of the phase field 

A superfluid flow with the velocity Vg — {h/m)'V^ 
constitutes a nonequilibrium excitation of the phase field 
<i>(r) in Eq. (|29|). At finite temperatures there will be also 
statistical excitations of the phase field. In the following 
we discuss such thermal excitations or fluctuations. This 
discussion leads to a qualitative argument for the exis- 
tence and the meaning of the coherence limit k^oh in Eq. 




We assume an order parameter field of the form 

*(r)- y^exp[ia>(r)] . (32) 

This means that we do not consider spatial variations of 
|\I'(r)| in this subsection. 

The average momentum of the fluctuations of the 
phase field is denoted by 

Ajfluct = (|V$|) . (33) 
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This momentum kfiud will be a function of the tempera- 
ture. 

For the single-particle states of the IBG we use real 
functions ip^. We may then consider the possibility of 
phase fluctuations for a low- lying state with rik 3> 1, i. e. 
phase fields $k(r) that are introduced by the replacement 



<Pk — > (pkexp[i<i>k(r)] 



(34) 



Let us first assume that these additional phases van- 
ish, $k ~ 0. In this case, the average kinetic energy 
h'^k^^^^^/2m of a condensed particle would exceed that 
of a noncondensed particle with k < fcfluct- The energy 
sequence of the single-particle states is, however, a pre- 
requisite of the BEG; the condensate must be formed by 
the particles with the lowest energy. In order to preserve 
the energy sequence of the low-lying states we require the 
phase ordering 



$k(r) = $(r) for k < fcfluct 



(35) 



This argument does not apply to the states with mo- 
menta above ftfiuct- 

It is plausible to assume the phase ordering Eq. ( |35| ) 
not only for statistical but also for nonstatistical excita- 
tions. This means that the coherence limit introduced in 
Eq. (^ should be identified with the average momentum 
of the fluctuations: 



(36) 



In a flrst approach, the coherence limit fccoh in Eq- (|^) 
can be considered as a parameter that has to be deter- 
mined from the fit to the experimental superfluid density. 
In contrast to this point of view, we present now theo- 
retical arguments for the leading asymptotic behavior of 

^coh ■ 

The critical exponents of gs and fcfluct may be derived 
from the following assumptions: 

1. The Landau part Gl of the thermodynamic poten- 
tial has its well-known form (compatible with the 
IBG). 

2. The fluctuation part Gsuct contains the extra fac- 
tor Qs/ go due to the assumed comotion of noncon- 
densed particles (as compared to a standard GL 
functional) . 

3. The fluctuation part scales in the same way as the 
Landau part. This scaling invariance is the decisive 
physical assumption. 

We write Gql = / (Pf 3gl for the thermodynamic po- 
tential. The first two assumptions lead to an integrand 
of the form 



ffGL = fffluct + gh 



2 m go 



In contrast to Eq. (g|), we do not introduce an as- 
sumption about the temperature dependence of the ratio 
gs/go- The present discussion will rather deduce this 
temperature dependence. 

Asymptotically, the Landau part behaves like 



(38) 



The third assumption, the scaling invariance, implies 
that ((^fluct) must have the same leading t dependence. 
Because of (|V^'p) = Pofcfluct have (gfluct) oc gsk^^^^. 
If this term scales with the critical exponent of kd^ct 
must be smaller than 1. This implies that gcoh oc /ccoh = 
^fluct is asymptotically large compared to po '^^^ 1^1- There- 
fore, 



Qs = go + Qcoh Pcoh ~ fcflu 



(39) 



The last step follows from Eqs. and (|36|). This leads 
to 



uct fluct 



(40) 



The scaling assumption (gnuct) ^ (sl) ^ implies then 

fcfluct ~ . (41) 

According to Eq. ( p9| ) the superfluid density has the same 
critical exponent, gs ^ fcfluct \t\'^^'^- 

To summarize: The assumption of the coherent como- 
tion (leading to gicoh oc fccoh = fcfluct) and the scaling of 
the fluctuation part ((gfluct) oc pcoh fcfluct 

oc fcg^^^J with 

the Landau part {(gh) ~ i^) determines the critical ex- 
ponent f — 1/3. 

IV. CORRELATION LENGTH 

A. Derivation 

We derive the correlation length in our effective GL 
model. For this purpose we add a standard coupling 
term — /i • 4" to the integrand of the energy functional 
( ^ ) where h is a. fictitious external field. The variation 
of the energy functional with respect to the ^' field leads 
to the field equation 



ao X'' 



27r/|i|V3 



A* - 2 



Acp 



'7 



+ h{r) = . 

(42) 



A small external field h{r) — dh\^exp{ik ■ r) will change 
the field from its equilibrium value (4') to ^E* = + 
(5\E'kexp(ik • r). We insert these expressions for h{r) and 
^'(r) into the field equation. In zeroth order in S'i' this 
yields the equilibrium value (I^Pp). In first order in 5"^! 
we obtain 



(37) 
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Acp ft 



27r/|t|i/3 



(5*k = Shk . 

(43) 



This susceptibility of the system with respect to an ex- 
ternal field is given by 



X{k) 



Shu 



x(o) 



1 + /fc2 ^2 



(44) 



The k dependence displayed in the last expression follows 
from Eq. (HSh. It implies 



(5*(r) cx exp(-r/^) . (45) 



h{r) = hoS{r) 



This means that the quantity ^ introduced in Eq. ([4^ ) 
defines the correlation length. Comparing Eqs. ( ^3[ ) and 
(H) yields 



27I-|t|l/3 



2^{t + 3 



Inserting the equilibrium values 

= f\t\/v for t < we obtain 



(46) 



for t < and 



m 



(i>0) 



where 



Co 



A. 



47r V Acp/fcE 



1.53 A 



(47) 



(48) 



The 



For the numerical value we used Eqs. (HJ) and (] 
result is independent of the parameter /. 

The * theory by GS leads to Eqs. (|4^) and (||), too, 
but with an additional factor [(3 + M)/3]^/^ for f+ and 
[(3 + M)/(3 + 3M)]i/2 for here M is an extra param- 
eter. For M = the results by GS and the present model 
coincide. For the favored M value by GS (M = 0.6 ± 0.3 
according to Eq. (pO|) in Ref. ^ the GS values differ from 
ours by about 10%7^ 

The validity of the result Eq. ( |47| ) is restricted to the 
temperature range \t\ < 0.01, Eq. ([2^). Outside this 
asymptotic region we expect the following behavior: For 
\t\ > 0.2 we find (Fig. [l]) ft ~ Po corresponding to i? « 1 
in Eq. (^). This implies that in the temperature range 
from \t\ = 0.01 to \t\ = 0.2 the correlation length will 
change its behavior from ^~ ~ |t|~2/3 ^- ^ |t|~i/2. 



B. Critical fluctuations 



our energy functional. The following discussion proceeds 
along the lines followed in a standard GL approach (see, 
for example, Ref. p^ . 

Equation ( ^3|) determines the deviation (5^* from the 
equilibrium due to an external field. Deviations from 
equilibrium occur also due to fluctuations. The fluctu- 
ations of the thermodynamic potential are of the size 
SG — ©(fceT). This determines the size of the thermal 
fluctuations i5^'thorm of the order parameter field. 

We write * = (^i) -|- S-^ and expand Ggl = (Ggl) + 
(5Ggl into powers of 6"^. The term linear in JvP van- 
ishes because we expand around the equilibrium value. 
Therefore, JGql starts with the quadratic term: 



SG. 



GL 



d^r 



1 



2a|«'|2 



(49) 



Here g denotes the integrand in Eq. (^3|) . The derivative 
has to be taken at the equilibrium value: 



= 2 



Acp 
fee 



2 Acp 



P 



(i>0) 
(i<0) 



(50) 



The typical range of a thermal fluctuation (J'I'thcrm is 
given by the correlation length £,(t). For the follow- 
ing estimate we may, therefore, replace the ^ cPr in Eq. 
( p9| ) by Using this approximation and (JGql, therm — 
0{kBTx), Eq. IM) becomes 



^Ggl e d^g 



kBTx 2 dm 



1(5* 



therm 



0{1) 



(51) 



Inserting Eq. (|5y) and omitting factors of the order 1 we 
obtain 



\S^ 



therm 



0(1) 

\tmr 



\t\ 



(52) 



For the last step we used Eq. f - |t|-2/3. The appli- 
cability of a GL approach is restricted to temperatures 
where the fluctuations are small (or at least not large) 
compared to the mean value of the field. The relevant 
ratio is 



the 



(1*1 



const. 



present model 
GL 



(53) 



In a standard GL approach the ratio diverges for \t\ 0; 
this excludes the application of the model in the asymp- 
totic region. In our effective GL model the ratio remains 
finite for \t\ 0. This is the same behavior as in a 
standard GL approach in d = 4 dimensions. 



A standard GL ansatz (in c? = 3 dimensions) breaks 
down for |t| ^0 because the ratio of the fluctuating 
to the equilibrium field diverges. In this section we will 
show that this problem does not exist in our effe ctive GL 
model. This is due to the scaling invariance (Sec. [II C) of 



V. TRANSITION TEMPERATURE IN HELIUM 
FILM 

We consider a helium film on the y-z-plane extending 
from X ^ to X = D. We look for a nonzero solution 
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^{x) of the field equation that obeys the boundary con- 
ditions 



*(0) = ^{D) = 



(54) 



This implies that a nonvanishing ^'(x) is inhomogeneous. 
Due to the kinetic energy term in Eq. ( |2^ ) such a non- 
vanishing 4'(a;) becomes possible only at a temperature 
T\{D) below Tx = T\{oo). We calculate this transition 
temperature T\{D) as a function of the film thickness D. 

In an experiment the plane a: = might be a solid wall, 
and the plane x = D may define the free surface of the 
helium film. In this situation one has a solid layer (thick- 
ness d) of helium at the wall so that the first boundary 
condition actually reads "^{d) = 0. This point amounts 
to replacing D hy D ~ d and is ignored in the following. 
For a discussion of the appropriate boundary condition 
at the free surface we refer to GS0. 

We consider the field equation Eq. ( ^ ) for h(r) = 
and t < 0. We introduce the quantity of Eq. (Es) and 
specialize to the one dimension of interest; 



dx"^ 



^ 1- 



= 



(55) 



For determining the onset of a nonvanishing ^P, the 
term may be neglected. The solution of the remaining 
differential equation is \E' = A sin(a; |tp/'^/^o + ct) where 
|Ap < f\t\/v- The boundary condition ^'(0) = implies 
a = 0, and ■^'{D) = leads to fltp^^/Co = nn where 
n is a positive integer. We are looking for the highest 
temperature T or the smallest \t\ value {t < 0) for which 
such a solution exists. That means we have to take the 
n = 1 solution for which the boundary condition yields 



-2/3 _ 



= D. 



(56) 



This defines the maximum T value for a nontrivial so- 
lution. For a more elaborate discussion of the under- 
lying mathematics we refer to Ref. 0. Using t{D) = 
[TxiD) ~ Ta(^)]/Ta(oo) we resolve Eq. H) for Tx{D): 



Tx{p)=Tx{^) 



1 - 



1 - 



D 



3/2 



4.81 A 
D 



3/2 



(57) 



Here Ta(oo) denotes the transition temperature of the 
bulk system. The numerical value follows from Eq. (48). 

As discussed by GS, the experimental and theoretical 
situation in a helium film is more delicate. According to 
GSu there is a caloric transition, and at a slightly lower 
temperature another transition to superfluidity. The con- 
dition for the second transition is somewhat different 
from that one discussed-in this section. As shown by 
Kosterlitz and Thoulessllj, superfluidity sets in when the 
superfluid density (averaged over the film thickness) ex- 
ceeds a certain minimum value. Evaluating this condition 



in our model leads to an expression of the form (p7|), too, 
but with a slightly different length parameter (instead of 
^Co). 

For the compaijison with the experimental findings we 
refer again to GSo. One finds indeed transition temper- 
ature shifts that are proportional to D^^^^, and length 
scales similar to the calculated one. 



VI. DENSITY PROFILES AT A BOUNDARY 
A. Introduction 

We consider a helium film with a thickness D much 
larger than the correlation length ^. At the wall the 
superfluid density will raise from zero to its bulk value 
Qs (bulk) (for simplicity we ignore again the solid layer at 
the wall). Similarly, the superfluid density will fall off 
again towards the free surface. 

The following discussion is restricted to the density 
profile at one boundary. This means that we determine 
how the density profile 



/(^) = 



gs(a^) 
gig (bulk) 



(58) 



raises from zero at a; = to its bulk value. This corre- 
sponds to the boundary conditions 



/(0) = and /(oo) = l 



(59) 



Third sound measurements determine the value of the 
superfiuid density averaged over the film thickness D. 
By = (bulk) (D — 2^hcai)/i^ one defines the healing 
length fhcaiCf) (for a more comprehensive discussion see 
§ 41 in Ref. 0). For a given profile f{x) this healing length 
may be calculated by 
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heal 



= lim [a— dxf{x) I 
\ Jo J 



(60) 



We approach the question of the density profile in two 
different ways: 

1. In the considered asymptotic regime {\t\ < 0.01) 
the density gicoh is the predominant contribution to 
the superfiuid density. The particles forming Qcoh 
have a distribution of momenta up to about fccoh- 
This leads to a healing length ^hcai i 1/ ^coh ^ 

|t|-2/3. 



2. We solve the field equation (55) with the bound- 
ary conditions Eq. (p9[). This yields the profile of 
condensate density go, Eq. (|30|), and to a heahng 

length ,ehcai,2~ Co i^r'/^. 

The composition (|^) of the superfiuid density leads thus 
to the prediction that the boundary profiles of the den- 
sities Qs and go are different. 
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B. Profile of the superfluid density 



D. Discussion 



For |i| <C 1 the superfluid density may be written as 
Qs ~ Qcoh oc / dkk^ (rifc) |¥'k(i")P ■ 



(61) 



The (/3k (r) are the real single particle functions of low- 
lying noncondensed states. For \t\ <^ 1 the occupa- 
tion numbers (^) are (rifc) « l/('«^ + t^). Because of 
^ k^^Yl may write k'^/{k? + t^) « 1. At the bound- 
ary plane a; = all single particle functions must vanish, 
i.e. |(y5k(i")P oc sin^(fca;a;). Since k"^ = + ky + k^ is 
limited by k^^^^, we assume that k^ is limited by k^^^^/3. 
Putting these (partly approximate) steps together yields 



Qsix) CX 



fccoh/\/3 



tkoh 

2 V3 



dkx sin^(fca;a;) 

sin {2kcohx/V3) 
Ax 



(62) 



Including all constants, the first term on the r. h. s. equals 
the bulk value; it is obtained as the limit of Eq. ( |62|) for 
large x. Therefore, we may read off the density profile: 



Qsix) 



1 



sm 



(2fccoha;/V3) 



/i(x) . (63) 



£is(bulk) 2kcohx/V3 
We evaluate the healing length ( |60| ) for this profile 
TT \/3 a/3 a, 



Cheal, 



4fceoh 2C(3/2)ao 



|ir2/3 =o.82A|t| 



2/3 



(64) 



We used Eqs 
Eqs. (M) and (|2 




[13| ) and (|lj) for the second step, and 
for the numerical value. 



C. Profile of tiie condensate density 

Using the dimensionless quantities = 'i'{x)/'ii(hu\k) 
with ^'(bulk) = ^/f\t\/v and y ^ (x/^q) \t\^^^ , the field 
equation (|5|) reads 

r{y) + Hy)[^-Hyf]=o ■ (65) 

The boundary conditions Eq. (p3) become "0(0) = and 
"0(00) = 1. The solution of Eq. p5| ) obeying these bound- 
ary conditions is 7/;(j/) = tanh(y/\/2). Returning to the 
original variables we obtain the profile of the condensate 
density (pO|), 

ga{x) 



£io(bulk) 



= tanh 



V2eoit|-^/^ 



(66) 



Inserting this profile into Eq. 
length 



J) yields the healing 



&eal,2-^Co|tr'/^ = 2.16A|t|- 



-2/3 



(67) 



This differs by a factor of 2 from the correlation length 

r(0 • 



The calculated density profiles of the superfluid and 
condensate density, Eqs. ( p^ ) and (|66|), respectively, are 
depicted in Fig. ^ These results refer — as most results 
in this paper — to the asymptotic region \t\ < 0.01, Eq. 




FIG. 2. Superfluid density Qs{x) and the condensate den- 
sity Qo{x) as functions of the distance x from a boundary. The 
distance x is measured in units of Co|i|~^''^. For large x the 
densities approach their bulk values. 



The occurrence of different profiles does not alter the 
calculation of Sec. 0: For a superfiuid flow we must have 
a nonvanishing order parameter field. The condition ob- 
tained in Sec. ^ is sufficient to accommodate both, qq 
and gcoh, within the film. 

The experimental result for the healing length cited 
in Eq. (41-47) of Ref. is Cai ~ 1.55 A |i|-2/3. The 
question is which one of the calculated results, (^) or 
(|67|), should be compared with the experimental one. 
Both results scale correctly. Relative to the experimen- 
tal prefactor, one theoretical prefactor is about 50% to 
large, in the other case about 50% too small. Because 
of gcoh(bulk) > eo(bulk) (for \t\ < 0.01) one might be 
inclined to consider the correlation length ( |6^ as the 
relevant quantity. However, (i) the macroscopic 5" field 
defines the macroscopic phase field necessary for superfiu- 
idity, and (ii) the simple averaging procedure ( |60| ) might 
not be appropriate for a profile like fi{x). In view of 
this more involved situation, we restrict ourselves to the 
statement that the theoretical results are in a reasonable 
vicinity of the experimental one. 

A direct experimental verification of the predicted pro- 
files shown in Fig. ^ appears hardly feasible. The pro- 
file may have some influence on third sound calculations 
the results of which might then be tested. Jjhe conden- 
sate density is already in the bulk systenJlJ a hard to 
access quantity. 

The different profiles calculated here are obviously in 
variance with the GS approach. It is a future task to find 
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differing predictions of our model that might be readily 
subject to an experimental test. 

VII. CONCLUDING REMARKS 

We have proposed an effective Ginzburg-Landau model 
that is compatible with the bulk properties of liquid he- 
lium. We have investigated the basic properties of this 
model and its most simple applications for finite geome- 
tries. 

Our model is based on the IBG and the assumption of 
a coherent comotion of condensed and noncondensed par- 
ticles. The occurrence of fractional powers in the energy 
functional can be made plausible on this basis. The foun- 
datioruci our model is thus different from the 5" theory 
by Gffltl. 

The results obtainecLso far are mostly similar to that of 
the 4" theory by GSETq. This is to be expected for other 
applications, too. The results for the density profiles at 
a boundary are, however, novel and in variance with the 

theory. 

There is a variety of further problems that might be 
treated in our approach (and that have already been 
treated in the ^ theory), like the influence of external 
fields or of impurities (ions) and the description of su- 
perfluid flow in various geometries. _ 

The ^ theory has been generalizeda in order to incor- 
porate (i) time-dependent problems in the field equation, 
(ii) dissipation effects and (iii) problems with a nonvan- 
ishing normal velocity. These generalizations correspond 
to additional terms in the functional of the thermody- 
namic potential. All these generalization can be intro- 
duced in our model in a straight-forward way, too. 
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